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ABSTRACT 

By  including  the  effect  of  fluid  loading  for  the  thin  spherical  shell  in  a  proper  manner,  so-called  shell 
theories  can  predict  the  water  borne  pseudo-Stoneley  waves  described  extensively  in  the  literature.  Shell 
theories  give  reasonably  good  results  for  the  motion  of  a  bounded  elastic  shell  by  using  the  assumption  that 
various  parts  of  the  shell  move  together  in  some  reasonable  manner.  Without  proper  fluid  loading,  however, 
shell  theories  do  not  predict  the  pseudo-Stoneley  resonances  observed  in  nature  and  predicted  by  exact  theory. 
With  proper  fluid  loading,  as  well  as  rotary  inertia  and  translational  and  rotary  kinetic  energy  terms,  a  shell 
theory  can  exactly  predict  these  water  borne  resonances.  These  resonances  are  predicted  by  the  shell  theory 
and  compared  with  results  from  exact  elastodynamical  calculations. 


1.  INTRODUCTION 

Interest  in  sound  scattering  from  submerged  elastic  shells  is  broadly  based.  For  example,  the  structural 
engineering  community  has  a  considerable  interest  in  the  area  of  nondestructive  testing  via  ultrasonic 
scattering,  while  in  the  area  of  remote  sensing  or  object  identification  scattering  from  objects  proves 
invaluable.  The  acoustical  scattering  from  elastic  objects  or  more  simply  the  generation  of  resonances  from 
bounded  elastic  shells  can  be  described  by  acoustic  and  elastodynamic  theory.  Exact  solutions  for  scattering 
from  elastic  targets  exist  for  target  shapes  for  which  the  elastodynamic  equation  is  separable.  For  three 
dimensional  targets  separability  is  only  possible  for  spherical  and  rectangular  targets.  For  more  complicated 
targets  numerical  techniques  must  be  used.  For  axisymmetric  smooth  targets  the  extended  boundary  method  of 
Waterman  has  proven  useful  though  limited  and  an  approximate  theory  would  be  desirable  for  more 
complicated  targets.  The  purpose  of  this  work  is  to  develop  an  approximate  theory  based  on  shell  theory  that 
would  be  of  use  for  general  shapes  that  could  aLso  include  structural  loading.  Shell  theories  afford  a  very 
powerful  methodology  to  "build  in"  structural  features  usually  via  a  variational  principle  in  which  some 
Lagrangian  or  Hamiltonian  is  constructed  by  introducing  physical  features.  The  simplest  theories  which 
always  assume  correlated  motion  of  the  inner  and  outer  surfaces  of  the  target  include  kinematic  features  and 
potential  energy  terms  based  on  the  generalized  Hooke's  law  which  leads  to  the  lowest  order  symmetric  mode. 
One  can  add  rotational  inertia  to  the  Hamiltonian  which  allows  an  antisymmetric  mode.  One  goes  on  from 
there  to  include  higher  order  corrections.  Further,  proper  fluid  loading  -  which  from  our  experience  is  no  ea-sy 
task  -  must  be  introduced  correctly.  The  proper  inclu.sion  of  fluid  loading  introduces  a  fluid  borne  wave:  the 
pseudo-Stoneley  wave.  Clearly  thin  shell  theories  are  rather  geometrical  and  must  be  constructed  for  each 
shape  due  to  their  dependence  on  shape  dependent  dynamic  factors.  Non-spherical  three  dimensional  object.'' 
present  a  problem  which  we  seek  to  address  in  the  future.  In  this  work  we  wish  to  show  the  development  and 
test  the  predictions  of  a  suitable  shell  theory  for  spherical  shells  that  can  be  readily  generalized  to  spheroid.'' 
and  finite  cylinders,  and  .specifically  we  wish  to  show  that  proper  fluid  loading  in  an  shell  theory  does  allo" 
the  prediction  of  water-borne  p.scudo-Stoneley  waves.  In  the  following  section  we  outline  the  derivation  ot  a 
shell  theory  with  many  of  the  featured  described.  The  theory  is  then  used  to  compare  with  the  exact  methi'd^ 
and  specifically  the  exact  theory  calculations  of  the  water-borne  wave  mode.  The  results  arc  then  diseiisn 
and  future  work  is  de.scribed. 
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2,  DERIVATION  OF  EQUATIONS  OF  MOTION 

In  spherical  shells  membrane  stresses  (proportional  to  p)  predominate  over  flexural  stresses  (proportional 
to  p^)  where 

We  differ  from  the  standard  derivation  for  the  sphere*  by  retaining  all  terms  of  order  /3‘in  both  the  kinetic  and 
potential  energy  parts  of  the  Lagrangian  and  by  considering  the  resonance  frequencies  for  the  fluid  loaded  case 
to  be  complex.  We  note  that  this  level  of  approximation  will  allow  us  to  include  the  effects  of  rotary  inertia  in 
our  shell  theory,  as  well  as  damping  by  fluid  loading.  The  parameter  itself  is  proportional  to  the  radius  of 
gyration  of  a  differential  element  of  the  shell  and  arises  from  integration  through  the  thickness  of  the  shell  in  a 
radial  direction.  We  will  use  an  implicit  harmonic  time  variation  of  the  form  exp(-/ty/).  We  begin  our 
derivation  by  considering  a  u,v,w  axis  system  on  the  middle  surface  of  a  spherical  shell  of  radius  a  (measured 
to  mid-shell)  with  thickness  h,  as  shown  in  Fig.  1. 

2.1.  Lagrangian  Variational  Analysis 


Our  Lagrangian,  L,  is 


L  =  T-V-¥W, 


where  T  is  the  kinetic  energy,  V  is  the  potential  energy,  and  W  is  the  work  due  to  the  pressure  at  the  surface. 
The  kinetic  energy  is  given  by 


T  =  edxded(p. 


where  the  surface  displacements  are  taken  to  be  linear: 


,,  X,.  X  dw 

«.  =  (!  +  -)« — 

a  a  do 


w,  =  w. 


motion  of  the  spherical  shell  is  axisymmetric  since  the  sound  field  is  torsionless.  Thus  there  is  no  motion 
'n  the  v-dircction.  Substitution  of  Eqs.  (4)  and  (5)  into  Eq.  (3)  yields,  after  integration  over  x  and  <(> , 

'r  r’'  •  /ir/  ,  2x-2  h^..dw  ,  h^.,dw.-,  Jr  ,  2^-2i,a  r/:\ 

T  =  ;rp  I  sin0[( - 7  +  —  +  /ifl  )m  -2( - -  +  — )«-—-(-( - — )(^— )'+( —  +  ha  )w  ]dd.  (6) 

0.  d'  ;)fi  12  dd  12 


in  terms  of  /?, 
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T  ^  +  -0.6p*  +6P‘)u^HlSP*  +(p^  +  l)w^]sinddd,  (7) 

•0  /Vn  /^n 


<90  "  "  de 


where  the  first  and  last  terms  in  square  brackets  in  Eq.  (7)  are  associated  with  linear  translational  kinetic 
energies  and  the  middle  two  terms  are  associated  with  rotational  kinetic  energies  of  an  element  of  the  shell. 

The  potential  energy  of  the  shell  is 

1  fit  fix  fA/2  - 

9  Jo  Jo  J_,,2(<^«»eefi  +  tT,,e^,)(ar  +  fl)*sin0<irJ0<f0.  (8) 


where  the  nonvanishing  components  of  the  strain  are 


I  f  du  \  X  ( du  d'w 


=-(COt0M+w)  +  -4-COt0fM--^]. 
a  a  \  do  J 


and  where  the  nonzero  stress  components  are 


^ea  ~  j  ^2 


“  j  _  y2  ^^89 )> 


where  E  is  Young's  modulus.  By  substitution  the  potential  energy  becomes 


1/  1  f''r2''fA'2  E  1  x.du  X  d^w  ,  xdw^  ,2 

V  =  -  - Z-- - r  [(!  +  — )~Tr - +  w  ]+{cot0[(l  +  -)M - l  +  w} 

2J0J0  v2  (;f  +  a)2(^  a  de  ade^  a  add 

+2  v{cot  0[(1  +  -)u  - +  w}[(l  +  sin  ddxdd, 

a  a  ad  a  oB  a  ad  j 


which  after  integration  is 


,,  nEh  ^  ^  /I 

V  = - X  {(w'  +  — r  +(w  +  MCOt0)‘  +2v(vv  +  — )(w  +  i<cot6 

1  _  Jo  Sd  dd 

o2rydu  d^w^2  2  ^2  ...  dw^,dii  w 

+P^[(-z - ttT cot  0(m-— )  +2vcot 0(m-— )(- - —p)]} sin 06/0. 

00  00'  00  00  00  00' 


Terms  in  the  potential  energy  proportional  to  /?'  arc  due  to  bending  stresses. 

And  finally,  the  work  done  by  the  pressure  of  the  surrounding  fluid  on  the  spherical  shell  is  given  > 
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w  -  2  «a-  f  p  ivsin  Odd, 

(15) 

where  p,  is  the  pressure  at  the  surface. 

2.2.  The  Lagrangian  density  and  its  equations  of  motion 

>-™  to  a  Ugrongian 

~  1  -  V'  +  2  v(w  +  ~)(w  +  «cot  ^) 

with  corresponding  differenUal  equations 

Q  _  <9 L  d  dh  d  dh 

du  dedu,~Jt^'  (17) 


0  =  4^ _ ^  d  dh  ^  d^  Bh  _  d^  dL 

^  dSdw,  dtdw,  dSdtdw^'^d^'^' 

O  =  (i.e;}=)[^  +  e„,d|-(v.oo.=  ti)a]-/J=0-/,=  co,d|^^Kl+v)+/i>(v  +  co.^d)J^ 

-  ^1(1.8/?'  +  6)3'  +  l)i^  _  (1  , 

and  '  * 

■'’*  ^  - Kl  +  v)(l  +  ^' )  +  )3'  cot'  9))  * 

+cot9((2-v  +  col'9))3'-(l+v))a 

'^’^'^^’“‘®^+^'<>+''+“‘’e)|^-/S'cot9(2-y+cot'9)^-2(l+vV 


(19) 
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These  differential  e(]uations  of  motion  (19)  and  (20)  have  solutions  of  the  form 

«('))= 

■=o  dt] 


H'(n)  =  Xvv/*,(n).  (22) 

11=0 

where  77  =  COS0  and  P/77)  are  the  Legendre  polynomials  of  the  first  kind  of  order  n.  When  the  differential 
equations  of  motion  (19)  and  (20)  are  expanded  in  terms  of  Eqs.  (21)  and  (22),  we  obtain  a  set  of  linear 

^uations  in  terms  of  U„  and  W,,  whose  determinant  must  vanish.  We  consider  two  cases:  with  and  without 
fluid  loading. 

2.3.  Vacuum  case 

'^e  simpler  case  is  that  when  the  spherical  shell  is  surrounded  by  a  vacuum  such  that  there  is  no  damping. 

In  this  case,  the  pressure  at  the  surface  vanishes;  =  0.  The  set  of  linear  equations  the  expansion  coefficients 
must  satisfy  are 

0  =  [Q^l  +  6^'  +  1.8)3^)-(l  +  /3^),c](/,+[ft^(3i8'  +  1.8)3^)-/3V-(l+v)]W„,  (23) 

and 

0  =  -Aj(Kr- 3)/3'  - 1.8/3"  + 1  +  v]f/,  +[Q'(1  +  + 1.8)3") -  2(1  +  (24) 

where  £l  =  (t)alc^,  xr  =  v  + A„  - 1,  and  A,  =n(7i  +  l).  In  order  for  Eqs.  (23)  and  (24)  to  be  satisfied 
simultaneously  with  a  non-trivial  solution  the  determinant  of  the  system  must  vanish: 

0  =  Q^(l  +  6)3^  + 1.8^" )(1  +  2)3^  +  1.8)3")+fi' {(3^' + 1.8)3" )Aj(v -  3))S' - 1.8)3"  + 1  +  V] 

-[2(1  +  V)  +  /3'  k-A,  ](1  +  6)3'  + 1.8)3" )  -  (1  +  ^' )  v(l  +  2^'  +  1.8)3" )} 

+(1  +  ^')k12(1  +  v)  +  /3' vAj -  AJ(7c-3))3'  - 1.8^"  + 1  +  v](^'x:+ 1  +  v).  (25) 

Since  there  are  no  damping  terms,  the  shell  vibrates,  in  theory,  forever.  Thus,  the  normalized  frequency  fl 

can  be  ^en  to  be  real.  Equation  (25)  is  quadratic  in  Q',  thus  we  expect  two  real  roots  to  (25)  and  thus  two 
modes  for  the  motion  of  the  shell.  They  are  the  symmetric  and  antisymmetric  modes. 

2.4.  Fluid  loaded  case 

For  the  fluid  loaded  case,  we  must  consider  a  modal  expansion  of  the  surface  pressure  in  terms  of  the 
specific  acoustic  impedance  z„.  In  its  most  general  form  this  is 


p{a,0,(t>)  =  ^m«^(cos0)cosm0. 


where 
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! 

i 


2  =  ipc~ - . 

h:m 

(27) 

The  specific  acoustic  impedance  z„ 

can  be  split  into  real  and  imaginary  parts: 

z,  = 

(28) 

where 

1  Km  1 

(29) 

and 

"  0)  /t'(ka) 

(30) 

For  the  case  we  are  considering  of  axisymmeiric  motion,  the  surface  pressure  is  given  by 

?.(«)  = -iz.V.(<;ose),  (31) 

/|s0 


or 

P.CS)  =  -<B=W','".«(cos9).  (32) 

ffsO 

Use  of  Eq.  (32)  in  our  set  of  differential  equations  of  motion  (19)  and  (20)  yields  the  following  set  of  linear 
equations  for  the  expansion  coefficients  in  the  case  of  a  fluid  loaded  spherical  shell: 

0  =  [fi'  (1  +  6^^  + 1. 8^" )  -  (1  +  /?^ )  +  [Q^  +  ISP*  )-p^K- (I +v)]W,,  (33) 


and 


0  =  -A„ [( V - 3))3'  - 1. 8/3 V 1  +  v]U„  +[a' {\  +  a  +  2p^ +l%p*)-2{\+ v)  +  Qiy-p^KX„]W„,  (34) 


'here 


and 


a  =  — 

P,h 


y  = 


£JL 


(35) 


(36) 


ihc  determinant  of  Eqs.  (33)  and  (34)  must  vanish.  However,  in  this  instance  the  value  of  Q  must  be 
to  be  complex;  the  resonances  have  a  width  that  depends  on  the  damping.  The  result  of  setting  this 
^■'^inant  to  zero  is 
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0  =  +  6/3'  + 1.8/3"  ){1  +  a  +  2/3'  + 1.8/3"  )+Q'/y(l  +  6/3'  + 1.8/3" ) 

+a' {(3/3' + 1.8/3" )AJ( AC- 3)/3'  - 1.8/3"  + 1  +  vj 
-(2(1  +  V)  +  /3'k-A  J(l  +  6/3'  +  1.8/3" )  -  (1  )AC(1  +  a  +  2/3'  +  1.8/3" )} 

+Q[-;>(1  +  /3' )  V]  +(1  +  ^' )  ac[2(1  +  V)  +  /3'  acA„  )  -  A  J(  v  -  3)/3'  - 1.8/3"  + 1  +  v'l(/3'  ac  + 1  +  V'). 

Equation  (37)  has  four  complex  roots.  From  work  with  an  exact  modal  solution  to  the  problem,  we  exp 
two  roots  to  be  associated  with  the  symmetric  and  anti.symmetric  modes  of  the  shell.  We  expect  the  other  i 
roots  to  be  associated  with  a  water-borne  pseudo-Stoneley  wave. 

3.  CONCLUSIONS 

The  next  step  is  to  plot  the  roots  of  Eq.s.  (25)  and  (37)  to  compare  the  resonances  predicted  by  the 
models  with  those  given  by  exact  modal  expansion  solutions.  By  suppressing  a  and  y.  the  model  associai 
with  Eq.  (37)  reverts  to  the  vacuum  case  model  a.ssociated  with  Eq.  (25).  Similarly  suppression  of  factors 
P  in  Eq.  (25)  will  result  in  a  reversion  to  a  previously  derived  solution  '.  In  a  previous  paper-  we  ranked  li 
three  different  models  according  to  their  degree  of  physicality  and  compared  their  re.sults  for  various  relati’ 
shell  thicknesses  against  each  other  and  against  the  exact  results  of  the  modal  expansion  model.  We  ab 
considered  the  limitations  of  each  of  the  models  including  the  exact  solution,  as  well  as  those  of  shell  mode 
in  general. 

These  models,  tluid  loaded,  vacuo  ca.se,  and  membrane,  were  successively  less  physically  sophisticate 
and  gave  successively  le.ss  good  comparison  with  the  exact  (modal  expansion)  results.  The  interested  reader 
referred  to  that  earlier  work  for  the  comparisons  of  the  models.  We  shall  content  ounselves  with  showing  th 
results  of  our  best  model ,  that  corresponding  to  the  roots  of  Eq.  (37).  In  Fig.  2  we  see  the  thin  spherical  stei 
shell  dilaialional  (symmetric)  and  flexural  (antisymmetric)  mode  resonances  calculated  by  this  model  whic 
includes  proper  fluid  loading.  Here  and  in  the  succeeding  figures  thick  means  h/a  =  0.1;  thin  means  h/a  ■ 
0.01.  The  shell  material  is  a  generic  steel  with  density  =  7.7  times  that  of  water,  shear  velocity  =  3.2- 
km/s,  and  longitudinal  velocity  v,  =  5.95  km/s.  The  surrounding  fluid  is  taken  to  be  water  with  density  p  = 
1000  kg/m^  and  sound  velocity c\„  =  1.4825  km/s.  The  exact  and  .shell  theory  calculations  agree  well  for  tht 
dilatational  (symmetric)  rc.sonanccs  and  exhibit  a  marked  improvement  over  the  results  of  our  membrane  [Eq 
(25)1  and  vacuum  models  [Eq.  (37)  with  a  and  y  suppressed)  for  the  first  .several  shell  theory  symmetric 
mode  resonances.  This  is  due  to  the  inclusion  of  both  rotary  inertia  and  fluid  loading  in  the  model.  The 
flexural  (antisymmetric)  mode  resonances  .show  the  appropriate  behavior  on  this  rather  limited  size  parameter 
scale.  Next,  in  Fig.  3  we  have  a  plot  of  thick  .spherical  steel  shell  dilatational  (symmetric)  and  flexural 
(antisymmetric)  mode  resonances  calculated  by  shell  theory  with  fluid  loading.  As  in  the  vacuum  case  a.s  well 
as  for  the  membrane  model,  the  dilatational  (.symmetric)  mode  rc.sonanccs  compare  well  for  exact  and  shell 
theory  methods.  Finally  in  Fig.  4  we  show  a  compari.son  of  modal  re.sonancc  predictions  for  a  1%  thick  steel 
shell  (thin  shell)  for  the  low  size  parameter  (=  ka,  where  k  is  wavenumber  and  a  the  radius  of  the  midshell) 
region,  just  where  one  would  expect  to  sec  pseudo-Stoneley  resonances,  if  the  model  is  capable  of  predicting 
them.  The  models  compared  arc  our  vacuum  model,  our  fluid  loaded  model,  and  an  exact  calculation  ol  the 
pseudo-Stoneley  resonances  in  this  region.  As  expected,  the  vacuum  model  mis.ses  the  water-borne  wave.s 
completely,  but  our  fluid  loaded  model  does  predict  them,  and  predicts  them  exactly. 

For  future  reference,  the  flexural  (anti.symmetric)  mode  re.sonances,  as  calculated  by  our  shell  theory  "iin 
fluid  loading,  do  not  have  the  correct  a.symplotic  limit  for  large  size  parameter,  although  they  do  exhibit  the 
correct  behavior  for  lower  values  of  ka..  The  proper  correction  fi)r  the  large  ka  a.symptotic  resonance  heha\i>'f 
could  be  most  easily  found  by  including  the  shear  stress  distortion  along  with  a  Timoshenko-Mindlin'*  shaiv 
factor.  Inclu.sion  of  the  shear  distortion  in  the  potential  energy  would  make  the  flexural  mt)de.s  asynirt‘’‘‘-' 
the  coincidence  velocity  and  the  .shape  factor  can  be  adjusted  so  that  they  asymptote  to  the  Rayleigh  velocib 
expected.  Work  is  now  under  way  to  do  just  ihi.s. 
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